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Abstract 

We propose a discrete non-abelian symmetry group which enforces maximal 
atmospheric neutrino mixing, while #13 = and the solar mixing angle 812 remains 
undetermined; without finetuning, 812 will be large but non-maximal. Our extension 
of the Standard Model has three right-handed neutrino singlets and implements 
the seesaw mechanism. Furthermore, we have an enlarged scalar sector with three 
Higgs doublets and two scalar gauge singlets; the latter have masses and vacuum 
expectation values of the order of the seesaw scale. Lepton mixing stems exclusively 
from the vr Majorana mass matrix, where non-diagonal elements are generated by 
the vacuum expectation values of the scalar singlets. The model predicts a neutrino 
mass spectrum with 7713 > m.2 > mi, and the effective Majorana mass of neutrinoless 
PP decay is equal to 11111712/1713. 
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1 Introduction 



The KamLAND experiment jT] has confirmed that the solar neutrino deficit — for a recent 
review see Ref. |2] — is explained by neutrino oscillations with matter effects jlj playing 
a decisive role. Whereas the solution of the solar neutrino problem, the large mixing angle 
MSW solution, displays a large but non- maximal mixing angle 8 12 [2], the atmospheric 
neutrino problem with mixing angle 623 requires sin 2 2# 23 > 0.92 at 90% CL, with a best-fit 
value sin 2 26*23 = 1 [El- 
There are many models and textures in the literature which attempt to explain large 
(not necessarily maximal) atmospheric neutrino mixing — for reviews see Ref. [E] . However, 
the closer the experimental lower bound on sin 2 2# 2 3 comes to 1, the more urgent it 
becomes to find a rationale for maximal atmospheric neutrino mixing. Arguments have 
been given [3|HJ|H] that in this case a non-abelian symmetry is required — for models along 
this line see, for instance, Refs. (TU1 HH E21 EE] • 

In this letter we propose a model which employs a discrete non-abelian symmetry 
group G with eight elements. The essential ingredients of this model are the seesaw 
mechanism J3j and a pair of scalar gauge singlets xi,2 in a two-dimensional irreducible 
representation of G. Neutrino mixing stems solely from the Majorana mass matrix Mr 
of the right-handed neutrino singlets vr. The vacuum expectation values (VEVs) of xi,2 
generate off-diagonal elements of Mr, and atmospheric mixing turns out maximal. It is 
essential for our model that the VEVs of the scalar singlets be at the (large) seesaw scale. 
The present model has some similarities with the model of Ref. [TOj, where part of the 
symmetry group is given by the non-abelian group 0(2) [To] . 

This paper is organized as follows. In Section |2] we introduce the multiplets of the 
model and in Section El we discuss the horizontal symmetry group G. The Yukawa sector, 
the VEVs of the scalar singlets, and the neutrino masses and mixings are presented in 
Sections El and respectively; in Section El we digress to explain how m M <C m T may 
be achieved in a technically natural way in the context of our model. After discussing 
some of the predictions of the model in Section |H1 we come to the conclusions in Section 

2 Particle content and symmetries 

Our model has three lepton families, including three right-handed neutrinos. We generi- 
cally denote e, /1, and r by a. Thus, we have three left-handed lepton doublets D a , three 
right-handed charged-lepton singlets and three right-handed neutrino singlets u aR . 

In the scalar sector, we employ three Higgs doublets <f>x, 2 , and 3 , together with two 
neutral singlets Xi an d X2- For the sake of simplicity we shall assume the singlets to be 
real. 

We introduce three symmetries of the Z 2 type: 



J> 2 r) flips the signs of all multiplets of the r family and l} 2 makes a transposition of the 



7L\ 




Dt, t r , u tR , X2 change sign; 

Dp <-> D T , /Ir <-> Tr, V^r <-> V t r, Xl *-> X2, 03 -» -035 

VeR-, v^r, v rR , 1; e R change sign. 



(1) 
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multiplets of the /i and r families. The auxiliary symmetry l}^ serves the purpose of 
allowing for m M 7^ m T while preserving the appropriate form of the neutrino mass matrix 
A4 U , as we shall see soon. The auxiliary Z 2 is equivalent to the one employed in the model 
of Ref. JU]> whose charged-lepton mass generation is identical with the present one. 



3 The group 

The symmetries zjj and do not commute, and together they generate a non-abelian 
group G with eight elements: 

set (G) = {e, g, h, gh, hg, ghg, hgh, ghgh = hghg} , (2) 

where e is the identity element, and g and h satisfy 

g 2 = h 2 = (gh) 4 = e. (3) 

The group G has five inequivalent irreducible representations: 

2: *->(o -1 S) ; (4) 

l CTCr / : g — > al, h — > er'l, with a, a 1 = ±. 

The 2 is the only faithful irreducible representation. It is obvious that (D^, D T ), (//r, tr), 
(iSpR, v t r), and (xi, X2) transform as 2, while 4> 3 is a 

The irreducible representations of the group G have a single non-trivial tensor product: 

2(8>2 = 1 ++ © l + _ © 1_ + © 1 . (5) 

From the 2 we read off that G can be thought of as a group of symmetry operations in 
the plane comprising rotations with multiples of 90° and reflections at the x and y axes 
and at the axes obtained by rotating the coordinate axes by an angle of 45°. This group 
is isomorphic to the crystallographic point group -D4 [16]. 



4 Yukawa couplings 



The Yukawa Lagrangian determined by the multiplets and symmetries defined in Sectional 
is given by 



Cy = - y\D e u eR + y 2 (d^u^r + D T v rR 

-y?,DeZR<\>\ - Va (d^/jr + D t tr) 02 - y 5 (py.HR - D T r R 
+ Y ^rC' 1 i^RXi + VtrX2) + H.c, 

where fa = ir 2 (p\. With VEVs (oU°|o) = vj /y/2 (j = 1,2,3), one obtains pi 
I2/3V1I and 



(6) 



m, 



V2 



V2 



\y4V2 - y 5 V3\ 



(7) 



Also, the neutrino Dirac mass matrix is 

Mr, = diag (a, b, b) , with a = y\v\ j V2 and b = y\v\ j V2 . (8) 

The smallness of the electron mass m e is correlated with the smallness of the neutrino 
masses, since m e is proportional to \vi\ while the neutrino masses are, through the seesaw 
mechanism, proportional to \vi\ . 



Note that v = y|fi| 2 + |t>2| 2 + |t>3| 2 — 246 GeV represents the electroweak scale. 
There is also a Majorana mass term 

1 r . 



Cm ~ 2 



M*u T eR C- l u eR + M'* (^rC-'^r + ^ r C-^ t r)] + H.c. (9) 



Since in our model the charged-lepton mass matrix and Mb are both diagonal, the 
lepton mixing matrix has to result exclusively from Mr. Indeed, non-zero VEVs of Xi,2 
will render Mr non-diagonal. Furthermore, examining the Lagrangians of Eqs. © and (JDj), 
we find accidental Z 2 symmetries 1$ an d ZjJ} , defined in analogy to ly , cf. Eq. (JTJ): 

Z { 2 } : D e , e R , VeR, Xh X2 change sign; 
: Dp, n R , v^r, xi change sign. 

The properties mentioned in this paragraph are reminiscent of the maximal-atmospheric- 
mixing model of Ref. JTO] . The Z2 symmetries of the present model correspond to the 
U(l)L a lepton-number symmetries of Ref. [ID) ; in the neutrino sector, the former symme- 
tries are spontaneously broken by the VEVs of Xi,2 leading to non-diagonal contributions 
to Mr, as we shall see in the next sections, whereas the lepton-number symmetries of 
Ref. [ID] are softly broken by the matrix elements of Mr. We further remark that, through 
a change in the basis used for the representation 2 of the group G, one can construct a 
model completely equivalent to the present one, where atmospheric mixing originates in 
the charged-lepton sector and solar mixing originates in the neutrino sector. 



5 The VEVs of the scalar gauge singlets 

We write these VEVs as 

(0 |xi| 0) = W COS7 and (0 |x 2 | 0) = Wsinj, (11) 

with W > 0. In the scalar potential V the ^-dependent terms compatible with the 
symmetries of Eq. (^J) are 

3 2 

v = ^(xl + xl) + {xl + xl^T.pj^j + ^ixl + xl) 

J'=l 

+ (x? - xi) (v4h + ri*4<h) + A' (xl -xt)\ (12) 

where \i, pj (j = 1,2,3), 77, A, and A' are constants. This potential is similar to the 
ones used in Ref. [17] . Note that a term (X1X2) 2 is implicitly contained in V. Thus, 7 is 
determined by the minimization of 

/ (7) = \'W A cos 2 2 7 + Re (r]v* 2 v 3 ) W 2 cos 2 7 . (13) 
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Provided A' > and the coefficient of cos 27 is smaller than the coefficient of cos 2 27 by, 
at least, a factor of 2, the minimum of / is at 

Re (77^3) MA s 

According to the seesaw mechanism ^3], we assume that W ~ |M|, \M'\ ^> v. Then, 
I cos 27 1 is very small. From now on we shall make the simplifying assumption that 

W 

(Oixi|0> = <0|x 2 |0>= (15) 

Corrections to this relation are of order v 2 /W 2 , i.e. they are suppressed by the square of 
the ratio of electroweak over seesaw scale. 

The potential in Eq. (JT2jl is invariant not only under the symmetries of Eq. ((TJ), 
but also under the accidental symmetries of Eq. (fTUj) . 



6 A possible mechanism for m /i /m r <C 1 

The muon and tau masses are given by Eq. (J2J). We find nothing wrong with the finetuning 
needed for obtaining <C m T , since an analogous finetuning is also needed in other 
models, in particular in the Standard Model. Yet, one could argue that in this case we 
do not simply have to choose one parameter to be small, rather we have to choose two 
products of unrelated quantities such that in m M those two products nearly cancel. 

In order to give a natural reason why m M <C m T , we may add to our model the following 
additional symmetry: 

T: //fl -»• -fi R , 02 <-> h- (16) 

This symmetry yields 

2/4 = -2/5- (17) 
In this way the finetuning is confined to the VEVs: 



V 2 - V 3 



V2 + V 3 



It 



One may show that a potential V which is invariant under T, in addition to the 
symmetries of Eq. (0), will in general allow for vacua which do not break T spontaneously, 
i.e. in which the VEVs of 0° an d 03 are identical. Indeed, if we write 

— = = U\ , —= = ue cosa, — = = ue^sma, (19) 

with Mi and u real and non- negative (u\ + u 2 = v 2 /2), and a belonging to the first 
quadrant without loss of generality, we find that there is a range of the constants in the 
scalar potential for which a — (3 and a = 7r/4, i.e. v 2 = v 3 . Thus, it is possible to obtain 

= by imposing the symmetry T and by not allowing it to be spontaneously broken. 

Now, in order to obtain ^ but small compared to m T , it is enough to introduce 
a soft breaking of T in the scalar potential. It is easy to see that the only possible term 
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breaking T softly — but keeping intact the other symmetries of the theory, in particular 
the ones in Eq. ((TJ) — is 

/isoft (402 - 4fo) , (20) 

in which /i so f t is a real coupling constant with dimensions of squared mass, which may 
naturally be assumed to be much smaller than v 2 . It is then easy to show that in the 
vacuum one gets cos la oc /i so ft/ u 2 <C 1, i.e. |t> 2 | 2 — |t> 3 | 2 <C |u 2 | 2 + |^ 3 | 2 , hence ^ <C 
m T . 

Thus, in the context of our model, and contrary to what happens in many other models, 
in particular in the Standard Model, it is possible to introduce a natural mechanism, 
operating through a softly broken discrete symmetry, for <C m T . 



7 The neutrino mass matrix 

With Eqs. ©, ©, and (|15jl it is straightforward to write down the Majorana mass matrix 
of the right-handed neutrino singlets: 

/ M M x M x \ , w 
M R = M x M' , with M x = . (21) 
V M x M' J ^ 2 

Using the seesaw formula [Hj 

M v = -M^M^Md (22) 

for the Majorana mass matrix of the light neutrinos, it follows readily from the 7L 2 
invariance of M D and M R — see Eqs. (JHJ) and (|2*T|). respectively — that M. v has the structure 

DDI 

/ x y y \ 

M v =\y z w . (23) 
\y w z J 

Maximal atmospheric neutrino mixing and # 13 = immediately follow from this structure 
of M. v — see, for instance, Ref. [18J — since that matrix has an eigenvector (0, 1, — 1) T . In 
our model, the form of M. v is a consequence of the non-abelian symmetry group G. The 
solar mixing angle is not fixed by the mass matrix of Eq. (|23|). but without finetuning it 
will be large and non-maximal |lUj . 

More specifically, the symmetric matrix of Eq. (|23J) is diagonalized by a unitary matrix 
V as 

V T M V V = diag (mi, m 2 ,m 3 ) , (24) 

where the light-neutrino masses m\, m 2 , and m 3 are real and non- negative by definition. 
The matrix V is decomposed as 

cos 9 sin 9 

V = e i& Ue if} , with U = | -sinfl/v 7 ^ cos^/v^ l/y/2 |, (25) 

- sin 9/V2 cos9/V2 -1/V2 
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where 9 = #12 is the solar mixing angle and a, (3 are diagonal phase matrices. The above- 
mentioned eigenvector of M. v appears as the third column of U. We assume, without loss 
of generality, that 9 belongs to the first quadrant while m<i > m\. The solar mass-squared 
difference is Awq = Am^ = m 2 — m 2 . The physical Majorana phases are 2 (/% — /3i) and 
2 ((3% — /?2); the other phases in V, notably the phases in the matrix a, are unphysical. 



8 Phenomenological consequences of the model 

Our model does not allow for a fully general matrix Mr, since it requires (M#) = 0, cf. 



Eq. (j21j) . A short calculation shows that this condition is rewritten as 



1 2i& 



1 / sin 2 ^ g2i/3i + cos 2 ^ e24/32 _ J_ 

2 I m! m 2 m 3 



0. 



(26) 



This relation allows one to express 

\(m)\ = \{M V )J = \m lC os 2 9e 2i ^ +m 2 sm 2 9e 2i(i2 

the effective mass probed in neutrinoless (3(3 decay, in terms of the neutrino masses: 

m\m2 



I Ml 



(27) 



'28) 



m 3 



Equation ()26|) also states that, in the complex plane, it is possible to draw a triangle with 
sides of length sin 2 #/mi, cos 2 ^/m2, and Therefore, the sum of any two of these 

numbers must be larger than, or equal to, the third one. One such inequality is 



sin 2 9 cos 2 9 1 
+ > — 

rrii ^2 ^3 



(29) 



Since we have chosen the convention m 2 > mi, it follows from Eq. (J29|) that m 3 > m\. 
Thus, our model requires a neutrino mass spectrum with m x < m 2 < m 3 . The remaining 
two triangle inequalities are summarized as 



1 

m 3 



> 



sin 2 9 cos 2 9 



rrii 



m 2 



(30) 



We can rewrite Eq. (|3Up in the form 



nil 



m? + Am 2 



> 



where we have identified Am 31 

A ™atm- 



atm 



nin 



sin 2 9 — cos 2 9 



nil 



(31) 



m\ with the atmospheric mass-squared difference 



Since 9 is large and Am 2 tm ^> Am 2 &1 it is clear that Eq. ()31|) always holds for nil 3> 



Am 2 tm . Thus, our model allows for a quasi-degenerate neutrino mass spectrum. As a 
matter of fact, Eq. (J31)) is satisfied either when nil is larger than a certain value of order 
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\J 'Am^tm, or when mi is in the vicinity of a certain value of order J Am© for which the 
right-hand side of that inequality vanishes. For the sake of concreteness, let us take the 
best-fit values Am 2 tm = 2.5 x 10" 3 eV 2 and Am 2 = 7.17 x 10" 5 eV 2 , tan 2 fl = 0.44 0. 
Then, we obtain from Eq. (|31|) that either 

3.16 x 1(T 3 eV < m 1 < 7.09 x 10~ 3 eV (32) 

or 

mi > 1.64 x 10~ 2 eV. (33) 

The value of |(m)| increases monotonically with increasing mi. For the lowest mi in the 
range of Eq. ()32|) one obtains |(m)| = 5.7 x 10 _4 eV. This value is too low to be within 
the range of sensitivity of the forthcoming experiments on neutrinoless /3f3 decay — for a 
review see Ref. [T2|. However, for the range of mi in Eq. (|33|) . we have the lower bound 
|(m)| > 5.8 x 10~ 3 eV, which is not far from future experimental sensitivity. 

In our model all the Yukawa- coupling matrices are diagonal. In the low-energy sector, 
where the effects of Xi,2 are given solely by their VEVs and their respective contributions 
to Mr, the only source of lepton-flavor violation is neutrino mixing. Denoting the scale 
of the right-handed neutrino masses, i.e. the seesaw scale, by ttir, it has been shown in 
Ref. [20] that in such a theory, in the limit m# — > oo and in the case of two or more Higgs 
doublets, the lepton-flavor-changing neutral-scalar vertices do not vanish, i.e., there is 
no decoupling in the Yukawa sector in that limit; the non- vanishing one- loop corrections 
are generated by charged-scalar exchange and by the mixing of the heavy neutrinos. On 
the other hand, all the lepton-flavor-changing processes which have either photon or Z 
exchange, or proceed via a box diagram, are completely negligible, since their amplitudes 
are suppressed by l/m 2 R , and thus decouple in the above limit. Processes which are 
mediated by neutral-scalar exchange and remain unsuppressed by inverse powers of m# are 

— > fi^^fi - , r ± — > fi ± e + e~, r ± — > e ± e + e~, and jj, — > e ± e + e~ . Their branching ratios 
have been estimated |H| to be of order Y s / (16ir 2 G F M$) , where Y is a typical Yukawa 
coupling, Gp is the Fermi constant, and M is a typical mass of the five "light" physical 
neutral scalars in our model. Using the reasonable values 16tt 2 GfMq ~ 10 and Y ~ 10 -2 , 
such branching ratios are of order 10~ 18 j2Hj- However, due to the proportionality to Y s , 
a slight increase in Y could easily bring the branching ratio for ^ — > e ± e + e~ close to the 
sensitivity range of future experiments. 

9 Conclusions 

In this letter we have discussed an extension of the Standard Model in the lepton sector 
which employs a horizontal non-abelian discrete symmetry group G. With respect to 
the Standard Model, the new multiplets are three right-handed neutrino singlets, which 
allow the use of the seesaw mechanism, two additional Higgs doublets, and two real scalar 
gauge singlets xi,2- The G-multiplets are such that all the Yukawa-coupling matrices 
of the Higgs doublets are diagonal. Before \i,2 acquire VEVs, the mass matrix Mr of 
the right-handed neutrino singlets is diagonal as well. Thus, non-trivial lepton mixing 
originates in (xi,2) 7^ 0- Maximal atmospheric mixing and the vanishing of the mixing 



8 



angle 6*13 are linked to (xi) = (^2)0 = W/y/2. We have shown that the scalar potential 
following from G-invariance allows such a relation among the VEVs of the scalar gauge 
singlets, apart from corrections suppressed by (v/W) 2 , where v ~ 246 GeV. If W is 
of the order of the seesaw scale, then these corrections are completely negligible. It is, 
consequently, also natural that the physical scalars associated with xi,2 have masses of 
similar order of magnitude. Due to the form of the neutrino mass matrix in Eq. ()23j) . the 
solar mixing angle is free, but large and non-maximal if we do no employ finetuning. 

Since there are two scales in the scalar potential of our model, the electroweak scale 
and rriR, the scale of the right-handed-neutrino Majorana masses, the model suffers from 
the naturalness problem associated with radiative corrections to the masses of the light 
scalars. Those masses receive radiative corrections of order m\. A high level of finetuning 
is required in order to eliminate those large radiative corrections through renormalization. 
A possible solution to this problem, without destroying the nice features of the model, 
might be supersymmetrization. 

Our model requires a neutrino mass spectrum with mi < m>2 < m.3. In addition, 
there is a lower bound on mi of about 3 x 10~ 3 eV, if we use the best-fit values for the 
solar mixing angle and for the mass-squared differences Am 2 tm and Am^. Neutrinoless 
(3f3 decay is not necessarily within the sensitivity range of future experiments, but it may 
well be, cf. Eq. (J3*3J) . Among lepton-flavor- violating processes, the same applies to the 
decay ^ — > e ± e + e~, while / u ± — > e ± 7 is strongly suppressed; lepton-flavor violation is 
just the same as in the model with soft lepton-flavor breaking of Refs. [TP} ED]. 
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